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We address the problem of a magnetic impurity in a two dimensional f/(l) spin liquid where the 
spinons have gap-less excitations near the Fermi-surface and are coupled to an emergent gap-less 
gauge field. Using a large N expansion we analyze the strong coupling behavior and obtain the 
Kondo temperature which was found to be the same as for a Fermi-liquid. In this approximation we 
also study the specific heat and the magnetic susceptibility of the impurity. These quantities present 
no deviations from the Fermi-liquid ones, consistent with the notion that the magnetic impurity is 
only sensitive to the local density of fermionic states. 
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I. INTRODUCTION 



A large number of theoretical proposals for the 
low-energy description of spin-liquid phases consider 
fractionalized fermionic degrees of freedom, the spinons, 
carrying spin 1/2 but no electric charge, coupled to an 
emergent U{1) gauge field, a gap- less photon-like mode. 
The spinons are gap-less having either nodal points^ or 
a Fermi-surface. The former case arises naturally in the 
slave particle approach to the t — J modeP but also in 
other physical contexts such as the half-filled Landau 
level"^ and the description of metals at a Pomeranchuk 
instability^. It presents non- Fermi liquid behavior due 
to the strong interactions between the spinons and the 
gauge field that lead to a spectrum with no well defined 
quasi-particles. This phase has a number of remarkable 
thermodynamical and transport propertied^, for example 
at low temperatures soft gauge modes contribute to the 
specific heat with a term proportional to T"^/^ Isl. 

Magnetic impurities embedded in a parent material 
provide an experimental probe to the bulk proper- 
ties and can help to discriminate between possible 
candidate phases. Moreover in order to be observed 
experimentally the system itself should be stable to 
a dilute density of such impurities. In a Fermi-liquid 
an antiferromagnetically coupled spin impurity leads 
to the well known Kondo efi^eciPl characterized by a 
cross over from the low temperature strong coupled 
regime, where the magnetic moment of the impurity is 
completely screened by the bulk quasi-particles, to the 
high temperature regime, where the impurity suscepti- 
bility follows a Curie- Weiss law. This cross-over occurs 
near the Kondo temperature which is an example of a 
dynamically generated energy scale. Since the under- 
standing of the Kondo effect the study of impurities in 
different bulk phases has attracted much attentio rP^, 
in particular for bosonic^*' and algebraic spin liquidd^^^^. 

The purpose of the present work is to study the 



behavior of a magnetic impurity embedded in a U{\) 
spin liquid with a Fermi-surface. Being a charge insu- 
lator this system still presents a Kondo like behavior 
since the spin degrees of freedom are free to screen 
the magnetic impurity at low energies. The article is 
organized as follows: in sec|ll]we describe the model and 
give some details of the \/N expansion (sec. II A| and 



II B I , the specific heat and the local spin susceptibilities 



are respectively computed in sec. |II C| and |IID[ Finally 
in sec |III| we conclude discussing the implications of our 
results. 



II. METHODS 

Starting from the t — J model in 2-D, the action de- 
scribing the spin-liquid phase with a spinon Fermi surface 
coupled to a compact U{1) gauge field can be obtained 
within the slave-boson formalisrrP or using a slave-rotor 
representation^'^, when fluctuations around the mean- 
field solution are considered. We assume that due to 
the presence of a large number of gap-less fermions the 
system is deconfined i.e. one can consider a non-compact 
U{1) gauge theorjff^. The partition function writes as a 
path integral over the spinon grassmanian fields fa=±i 
and the bosonic gauge fields = (ao,a) with action 



SsL = / d^xY,{fH^r+^ao)f. (1) 
•' <j=± 

+ [(a - za) /t] . [{d + la) /,] - izaoj , 

where h is the microscopic lattice volume and m the 
spinon mass. The integration over the temporal com- 
ponent of the gauge field oq acts as an on-site chemical 
potential for the spinons enforcing fXfa ~ 1- We 
use the notation J' (Px = Jq dr J d^x. 

At X = the interaction with the magnetic impu- 
rity is given by Sk = Sserry + Jxb J dr S/(0).S, where 



2 



Sserry IS the action of the free impurity spin. Jk is the 
Kondo coupHng and S/ (0) = /j, (t, 0) ctq,^ (r, 0). 
Using a fermionic representation for the impurity 
spin S = cJ^cTfj^^Ca, this term writes expHcitly 
Sk = J dr {J^a^iidr +iX)ca - iX- JKbSf^o.S} , 
where A is an integration parameter inserted in order to 
enforce the constrain 4 '^a = 1- 



A. Large N expansion 

Perturbative expansions for the Kondo problem are 
plagued with infrared logarithmic divergences signaling 
the fact that, for low energy, the system flows to a strong 
coupled fixed point where the impurity forms a singlet 
with the bulk electrons. Even if resummation of the 
divergent terms is possible this method is not well suited 
to describe the low temperature phase. Alternatively the 
large N expansion reproduces the essential features of 
the Kondo effect in the strong coupling regime. However 
for temperatures of the order of the Kondo temperature 
Tx, where a cross over to the asymptotic free regime 
is expected, this technique becomes unreliable due to 
the violation of the occupancy constrain and instead 
predicts a continuous phase transitiorP^. Therefore our 
results are restricted to the low energy regime. For the 
U{1) spin liquid the large N expansion corresponds to 
the random-phase-approximation (RPA) used to obtain 
most of the physical predictions for this phasd^. Re- 
cently the validity of this method applied to this specific 
problem was questioned^ since all planar diagrams 
where shown to contribute to leading order. A possible 
resolution was proposed in^^ using a double expansion 
to control higher loop contributions and essentially 
recovering the RPA result. 

In order to perform a saddle-point expansion we 
generalize the above action to su{N) following the stan- 
dard procedurd^E^ESl the Pauli matrices a — {cti, 0-3} 
are replaced by the generators of su{N) t — {..., r", ...} 
with the index a = 1, .., N"^ — 1 and the coupling constant 
is rescaled Jk — > 7^ • The representation of the impurity 
spin is taken to be conjugate to the spinous one. Using 
the Fierz-like identity® the Kondo term writes 

Sk = ydr|^4(^r + iA)cCT-iAQ/ 

+ ^fe(/l(0)c.4//3(0))+JA'}, (2) 

and the last term of Eq. (nj) is now multiplied by Q/ 
defined such that bY,^ fl (x) (x) ^ J2a = Qf ■ 



FoUowinglSl, the interaction term is decoupled in- 
serting a bosonic Hubbard-Stratonovich field x = ^e*"^. 
The integration over cj) can be absorbed by a shift 
in A leaving a single real dynamical variable n. The 
integration over the fermionic degrees of freedom can 



then be performed and the partition function writes 
Z = / DaDXDne-^ where 

= -^Trln[-G-i]-lTrln[-^^-i] + 



+ I dxo\ -tXQf\- I d^'x ^iao, (3) 



1 

J^' 



is the action for the bosonic fields a^, A and k only and 
and are the inverse of the full interacting prop- 
agators of the impurity and spinon fermions. We proceed 
performing a saddle-point expansion in the large N limit 
imposing the a static ansatz 



Mr) 



Kq ^ 0, 
6a, f_, ifJ- 



(4) 
(5) 
(6) 



At T = 0, the variations of the action in order to 01^2 are 
trivially zero and the ones for oq, A and k give, respec- 
tively. 
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(7) 
(8) 

(9) 



The first equation fixes the chemical potential ^, 
where G'o(iw„,k) = (ia;„ — ek) ^ is the bare prop- 
agator of the spinous with single-particle energies 
£k — 2?^^^ ^ '^(0) = 1^ is the spinon density of 
states at the Fermi level and A is a high-energy cutoff 
for the dispersion relation. The last two equations 
were obtained by Read and Newns for the Coqblin- 
Schriffer HamiltoniafiiSj. In the limit where A is much 
smaller than the Fermi energy but much larger than 
the other energy scales the propagator of the impurity 
fermions is given by Fq (ia;„) = _^ +iAsgncj (where 
A — 7rn(0)KQ6) corresponding to a Lorentzian density of 
states p{v) = \ (^_^ (see Fig. |4|(a)). The saddle- 

point values £c and A are thus the resonance position 
and the hybridization width respectively. Identifying the 
phase shift of an bulk spinon scattered by the impurity 

(5/ (w) — tan^^ ( e ^Lj ) ' '81 is a particular example 

of the Friedel sum rule. Finally Eq. (|9| defines the 
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Kondo energy scale ksTK = y'e'^ + A'^ = Ae b^mjn ^ 

At zero order in 1/iV there is no influence of the gauge 

fleld in the dynamics of the impurity. 

A comment about the procedure is in order at this 
point. One could imagine starting with the bulk theory 
flxed point obtained in ref.'^, this would correspond to 
flrst renormalize the bulk system propagators and then 
introduce the impurity. However since l/N is the small 
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Figure 1: Propagators, vertices and external lines used to 
obtain the bubble-like diagrams. 



paxameter of our expansion entering in both the spinon 
and the impurity Hamiltonians it is natural to start with 
the bare bulk action. The equivalence of both results 
can be checked replacing the bare spinon propagator by 
the integrating one. 
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Figure 2: Bubble diagrams obtained summing over the 
fermionic degrees of freedom. Due to parity mixed diagrams 
and are evaluated to be zero as well as the one labeled hyZ. 
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B. Fluctuations 

Fluctuations due to the bosonic fields are obtained 
summing the fermionic bubbles in the RPA approxi- 
mation. Without the Kondo term {Jk — 0) the prop- 
agator D^i, = n^J of the longitudinal and transverse 
components of the Gauge field is given by the density- 
density and current-current response functions. Using 
the Coulomb gauge V.a = the longitudinal part is 
fully gaped IIoo — ^ yielding to screening, by the 
spinous, of a U{1) test-charge. Therefore one can safely 
ignore the dynamics of uq. The transverse component 

Hi j = (^Sij — '^^^ n is gap-less and results from the 

Landau damping of the collective transverse modes by 
the gap-less spinous. For Ifl^l < f_F(j we can write 



n(ifi„,q) 



I a, 



XQ 



(10) 



127rm 



where 7 = ^ and x ~ 

Using the diagrammatic rules of Fig. [T] the bubble-like 
diagrams, including transverse gauge as well as k and A 
fiuctuations, are given in Fig. |2]and are divided in pure 
impurity diagrams, mixed diagrams and gauge diagrams. 
The transverse component of the gauge vertex is such 
that q X ji = — ^q x k. 

The impurity diagrams corresponding to the fiuc- 
tuations of K and A were obtained in ref.l^ and are given 
explicitly in the Appendix |A 2[ 



Figure 3: (a) Pieces of diagrams giving zero by parity, (b) X 
Diagram contributing to the free energy. 

It is easy to prove that due to parity considerations 
all the diagrams including the pieces of Figj3]-(a) are zero 



(see Appendix All. These include the mixed diagrams as 
well as the ones labeled by Z in Fig. |2] This implies that 
at this order in the 1/N expansion the gauge propagator 
decouples from the impurity degrees of freedom and the 
infiuence of the impurity scattering enters only true the 
X diagrams contribution: 
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(11) 



bnlFa (iw„) Go («w„, k - q) Go (iw„, k - q') x 



1 (k X q)^ (k X q')j- 

\ 

m ra^ |q| |q'| 



[Go (jw„_„, k) + Go (iw„+„, k)] 



N {Simp + Sa), whcrC 



Taking into account all non-zero terms the bosonic action 
3|, developed at Gaussian order, writes now Sf, = so + 



(12) 



so = -Trln [-G^^] - Tr In [^F^^] 



1 

Jk 



Ql 

N 



^3 + |i-^Ji 



is the value of the action at the saddle-point, 
eludes the fiuctuations of the impurity degrees of free- 
dom (given the original Read and Newns paper and in 
Appendix A 2 ) and Sa is the action for the transverse 



component of the gauge field 
1 f dqdq' 



(2^) 



a(ifi„,q)a(ir2t,,q') x 
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[5 (q - q') n {iQu, <i)+X {in^, -q, q') 



where a — qj_a and Xij — (q_L)j ( q'j^ ) X 



C. Specific Heat Capacity 

We compute the specific heat considering the temper- 
ature dependency of the free energy 
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So term gives 
energy of the 

N'^Vn{0)T 

impurity term (c^'"'^' ^ N^p{0)Ty The 1/iV terms 

carry the free energy contributions from the bosonic 
degrees of freedom. For low temperature all internal 
(fermionic) propagators of Figj2] can be computed at 
T — and the temperature dependency is given by 
the bosonic degrees of freedonpSl. The first next to 
leading order correction due to the impurity bosons 
(proportional to TrlnF in Eq. 13 1 has been shown to 



give a correction to the impurity contribution to the 
specific heai!^. Defining Cv 



(imp) 



7i 



imp 



T one obtains 

7imp = — 1) ^p{0) which can be interpreted as the 
suppression of one of the N impurity degrees of freedom 
due to the existence of a constrain. 

Since the fiuctuations of the gauge and impurity 
factorize new phenomena can only arise from the 
X corrections to the propagator of the gauge. In 
a system with a dilute number of magnetic scatters 
this term would be of the order of the density of 
impurities, in this case of a single impurity it is simply 
proportional to l/V. It is therefore natural to expand 
Trln[n + X] = Trln[n] + Tr [H-^X] + .... The first 
term in the expansion is responsible for the gauge field 

, 2/3 



(gauge) 



CX V 



contribution to the specific heat Ci 

The correction to the free energy AF = ^^Tr [n~^X] 
is given by the diagram of Fig. |3]-(b). It is easy to prove 
that such contribution vanishes remarking that one can 
rewrite it as 



(14) 



2bKlFo (iuJn) Go (ia;„, k) Go (ia;„, k) E/ {iuj„) = 

2 

where (jcj„) oc — isgn (a;„) |w„| ^ is the spinon self- 
energy^^ given in Fig. H-(b). The vanishing of such term 
is a consequence of the independence of from the 
spinon momentum. 



Thus the only correction to the specific heat due 
to the presence of the impurity is given by a correction 
to ')'(™P) since all other terms vanish either by parity 
considerations of by the above argument. 



D. Spin Susceptibility 

In this section we consider the local spin-spin 

correlations at the impurity site and its dif- 

coming from the impurity- 

= (S(r).S(O)) , impurity- 



ferent contributions 
impurity Ximp,imp (t) 

spinon Ximp, spinon 



from the 

(S/,(x=0)(t) 



local 



(S(r).S^,(,^o)(0)) and 



spmon-spmon X\ 



spmon.spmon 



■S/,(x=o) (0)) susceptibilities. In order to 



investigate the role of the impurity and gauge degrees 
of freedom we consider the 1/A^ corrections of the 
propagators and external vertices. 



(a) 



(b) 




Figure 4: Dyson's equation for the impurity (a) and spinon 
(b) propagators up to order 1/N, for sake of clarity symmetry 
related diagrams are not shown. 

Fig|4] shows diagrammatically the impurity and spinon 
propagators up to order 1/A^. One can see that to this 
order the impurity propagator has no corrections due 
to the presence of the gauge field since terms like 
f dk GoGo^f = vanish as a consequence of the 
independence of from the spinon momentum. Alter- 
natively one can use the renormalized spinon propagator 
to compute the self energy of the impurity (second 
term of Fq^ in Fig|4]-(a)) which would correspond to 
a rearrangement of the terms in |^(a) leading to the 
same result. The impurity propagator is thus the same 
as if the bulk was a regular Fermi-liquid. In this case 
one can use the results in ref.^^ where the fiuctuations 
of the bosonic impurity fields A and k, were shown to 
renormalize the Kondo temperature. 
Besides the self energy term the spinon propagator, 
given in Fig. |4]-(b), has also a 1/V^ contributions from 
impurity scattering, these can however be safely ignored 
in the computation of the local susceptibility since it 
would give a correction. 
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l/V" corrections due to 
impurity scattering 



Figure 5: Impurity-impurity and local spinon-spinon spin 
susceptibilities. The vertex corrections are given to order 
1/N , the bosonic propagators are obtained inverting the bub- 
ble like diagrams of Fig. [2] The external zigzag line carry 
spin and frequency indices. 



III. DISCUSSION 



We considered the Kondo screening in a bulk system 
of spinous strongly interacting with a J7(l) gauge field. 
While it is remarkable that Kondo screening can occur 
for a charge insulator, the results obtained here predict 
that no particular signature due to the presence of the 
gauge field can be measured if only the impurity degrees 
of freedom or local magnetic properties are monitored 
at the impurity site. 

The presence of the impurity destabilizes spin-liquid 
phase locally and Friedel-like oscillations are expected 
once the density of spinous is locally disturbed. This 
is due to the last term in |4]-(b), however they would 
equally be present if the density of spinous was changed 
by non-magnetic means as for example at the sample 
edges or near non-magnetic impurities. Such oscillations 
should be enhairced by the presence of the gauge 
fields however they are non-local measures. Local 
probes will be incapable of distinguish the bulk system 
from a Fermi-liquid. In particular the Wilson ratio 
R — ^"^'p,""vi^) _ _N £ ^]^jg ^g^gg jg ^YiQ same as for 

J(J+l)7i„p N-1 

a magnetic impurity embedded in a Fermi-LiquicP^. 



The impurity-impurity susceptibility is given at 
leading order in by the bubble diagram of Fig[5]-(a) 
(first term in the r.h.s.). vertex correctioirs due 

to the gauge field arising in the impurity-impurity 
susceptibility also vairish (see Figj5| since they coir- 
tain the terms like the ones in Fig|3]-(a). One thus 
concludes that the impurity-impurity susceptibility 
Ximp,imp has uo Contribution from the gauge field at 
this order in \/N. So the impurity degrees of freedom 
see only the local density of the spinous, in particular 
the result given irP^ for the static susceptibility hold: 
Ximp,imp («f^n = 0) = \j[J + 1)(2J -f- l)p (0) where 
TV = 2J+ 1. 

Gauge contributions are known to enhan ce Friedel- 
like oscillations in U{1) spin-liquid d ^ 1 1 1 '' l ^ .fll this is a 
consequence of the renormalization of the 2k p compo- 
irent of the susceptibility vertex. One could thus expect 
that the local spinon-spinoir susceptibility carried some 
trace of this behavior. Remarkably no vertex corrections 
to the local susceptibility due to the gauge field are 
possible since simple parity arguments like the one used 
in Appendix |A 1| show that the contribution given by 
the second diagram in the r.h.s of Fig[3]-(b) vanishes. 

Fiirally the crossed impurity-spinon susceptibility 
can also be shown to remain uiraffected by the presence 
of the gauge field using the same simple arguments. 
This shows that the local measurements of the suscep- 
tibility at the impurity site are insensitive to the gauge 
degrees of freedom. 
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Appendix A: Some details 
1. Prove that the diagram of Figj3]-(a) is zero 

The diagram of Fig. |3]-(a) is given by 

/dk (k X q) 
- — -2 Go (iWn, k) 1— r Go {iuin-u: k - q) . 
(27r) m\q\ 

Changing variables k — > k' = 2^q — k, where k' is 
obtained reflecting k on axes q (see Figj6|, leaves the 
norms |k| = |k'| and |k— q| = |k' — q| invariant and 
chairges the sign of k x q = — k' x q. Since Go («a;„, k) — 
Gq {iuJn, |k|) for a spherically symmetric Fermi surface it 
follows that W = -W = 0. 



6 



/ 




Figure 6: Change of variables that implies that the diagram 
in Fig[3]-(a) is zero. 



2. Impurity Fluctuation 

The impurity action at Gaussian level is given by 



->imp 



-y- 



A(ifi„) 



where 



O S O S 

SrSk SXSk 
SXSk SXSX 



is the fluctuation matrix. If one evaluate the fermionic 
Matsubara sums at zero temperature its entries are given 

by 



SxSkS 
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